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Abstract 

Using lattice approximations of M.'^, we develop a way to approxi- 
mate stable processes that are represented by stochastic integrals over 
K'^. Via a stable version of the Lindeberg-Feller Theorem we show that 
the approximations weakly converge as the mesh-size goes to zero. As 
an application, we improve upon previous approximation schemes for 
integrals with respect to linear fractional stable motions. 
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1 Introduction 

Stable integration is an important tool in the theory of a-stable pro- 
cesses. Similar to the theory for Gaussian processes, it is known (^ jST941 
Sec. 13.2]) that all stable processes, satisfying mild conditions, can be 
constructed from integrals of the form 



where is an independently scattered a-stable random measure on 
the measurable space {E,£) with control measure m and {ft)teT is 
a kernel such that ft G L°'{E,£,m) for a\\ t e T. If T = M and 
/* = l[o,t] then Xt is an a-stable Levy motion having independent 
and stationary increments (the symmetric case is the stable analog of 
Brownian motion). 

In this work, we approximate the finite-dimensional distributions of 
([T]) using a Riemann sum- type scheme. These approximations are use- 
ful for the dual purposes of intuition and simulation of stable processes. 
The weak convergence of our scheme is facilitated by a Lindeberg-Feller 
type stable limit theorem, which we have not previously seen in the 
literature. 

A couple of different discrete approximations of stable processes 
have appeared previously in the literature. One approach is Lepage's 
series which was improved upon in a series of papers by J. Rosinski 
(see |Ros01| and the references therein). In the present paper, we 
use a lattice approximation of stable integrals which extends, to / € 
L°'{M.'^), the "moving-average" discrete approximations of L-FSMs in 
[Dav70[ iMaeMl I Ast831 IDR85] corresponding to the case /* = l[o,t]. The 
work of |KM88| improved upon these earlier papers to obtain discrete 
approximations of slightly more general stable processes, while |AT92] 
showed that tightness of discretized L-FSMs cannot be achieved in 
the Ji-Skorokhod topology. In }KT95j . it was shown that discretized 
L-FSMs satisfy the fractional ARIMA equations and a closer look at 
issues concerning absolute convergence was taken. 

A secondary purpose of this work is to generalize certain Gaussian 
integrals to the a-stable case and, as in jKM88| . we then approximate 
such integrals with the scheme just described. In the past fifteen years 
or so, there has been an effort to develop stochastic integrals with 
respect to a broader class of Gaussian processes than just Brownian 
motion. In particular, consider Gaussian processes with stationary 




t e T, 
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increments, but replace the independent increments condition with the 
weaker condition of self-similarity. Normalizing the variance at t = 1 
to unity, one gets the single parameter family of fractional Brownian 
motions (FBM) with Hurst self-similarity parameter < H < 1. 

The theory of integration with respect to FBM is difhcult because 
FBM is not a semi-martingale. Nevertheless, rapid progress has been 
made using several different approaches (with significant overlap be- 
tween them). Roughly speaking, they can be categorized into four ap- 
proaches which use, respectively, fractional derivatives and integrals, 
Malliavin calculus, fractional white noise theory, and path- wise inte- 
gration (see IBH0ZO81 ). 

In Section[3l we consider a generalization of the FBM integral based 
on fractional integro-differentiation to a-stable analogs of FBM called 
the linear fractional stable motion£0 (L-FSMs). By "a-stable analog", 
we mean that a L-FSM is a self-similar, symmetric stable process with 
stationary increments. Any process with these properties is called a 
FSM. In contrast to the Gaussian picture, for each admissible {a, H) 
pair, there is not a unique (normalized) FSM, up to finite-dimensional 
distributions. Moreover, for each (a, H) pair with < i? < 1 and 
H ^ 1/a, there are infinitely many L-FSMs. These L-FSMs are rep- 
resented by (U) where _E = R is equipped with Lebesgue measure. Ma 
is symmetric, and 

/r'(^) := (2) 
a ((t - xf^-'/-^ {-x)1-''") +b({t- xf_-"'^ {-xf_-"'^) 

for properly normalized order pairs (a, b) where a, 6 > (see |ST941 Sec. 
7.4] for more details). Here x_ = |a;| if a; < and otherwise (similarly 
for a;_|_). The family of L-FSMs were the first FSMs to be constructed 
and studied, and much is known about them. Our motivation comes 
partly from |PTOO| which handles the a = 2 case. As in their work, we 
restrict ourselves to deterministic integrands, but |PTOOI shows that 
even in the a = 2 case, the theory for deterministic integrands is not 
completely trivial. 

An integral with respect to L-FSM will be defined as an integral 
with respect to a linear fractional stable random measure which we 
define for a > 1 and all permissable Hurst parameters < H < 1 . We 
have recently learned that when H > 1/a, jMS08| has developed sim- 
ilar integrals and also discrete approximations for them. However, the 
convergence results for their approximations concern a strictly smaller 
class of integrands. In particular, they require bounded integrands 
which are piece- wise continuous (we require no continuity or bounded- 
ness) and which must satisfy a faster tail decay than ours. 

^The term linear fractional stable motion was introduced in [CM89] due to its close 
relation to linear time series (moving average processes). 
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The rest of the paper is organized as follows. In Section 2, we 
review the notion of stable random measures and present our result 
concerning the convergence of discretizations of stable random mea- 
sures. In Section 3, integrals with respect to L-FSMs are defined, and 
their approximation by moving averages of i.i.d. random variables are 
discussed. Section 4 is devoted to the proofs. 

2 Discrete approximations of SaS random 
measures 

A useful viewpoint is that a random measure is a stochastic process: 

Definition 1 (Random measure). Let {E,£) be a measurable space 
and V be a vector space of measurable functions / : i? — >■ R. A ran- 
dom measure on {E,£) is a stochastic process (M[/])/gy satisfying the 
linearity property: for all a^^a-i G M and fi, f2 G V , 

M[aifi + a2/2] = aiM[fi] + a2M[f2] almost surely. (3) 

Let us make a few comments concerning this definition. First of all, 
the linearity property ([3]) ensures that the finite-dimensional distribu- 
tions of the process (M[/])ygy are determined by its one-dimensional 
distributions. E 1a ^ V for A e S, we note M{A) = Af[l^] which is 
thought of as the random measure of the set A. If M (Ai) are indepen- 
dent for disjoint sets Ai,. . . ,Ak, then M is said to be independently 
scattered. For general f £ V, to emphasize the analogy with usual in- 
tegration, the notation M[f] — f{x)M{dx) is often used. Finally, if 
one so pleases, one may also view the random measure M as a random 
linear functional on the linear space V (see for example |Dud69| ). 

Let <S„(cr) be the symmetric a-stable (SaS ) law of index a G (0, 2] 
with cr > being the scale parameteiH. We denote the characteristic 
function of Sa (c) by 

A„(6i) =exp(-|CT6i|"), 6'gR. (4) 

To reduce notation, when cr = 1 we simply write Sa — 5^(1). 

We now consider the class of independently scattered SaS random 
measures, i.e. those where M[f] is SaS for all f £ V. Suppose that 
{E, £, m) is a measure space where m is a cr-finite measure and £q is the 
class of measurable sets with finite m-measure. Following |ST941 Sec. 
3.3], we say that the independently scattered SaS random measure Mq, 
has control measure m if Ma{A) has distribution Sa(m{AY/°') for all 

^Our approximation in Thm l2.l1 as well as the Lindeberg-Feller result, can be extended 
to stable distributions with skewness u ^ Q, however, to simplify calculations and notation 
we have assumed symmetry. 
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A G Eq. For such random measures, it can be shown that V — L"{E) 
(see |ST94[ Ch. 3]) and that the distributions Ma{A),A e Sq uniquely 
determine the characteristic functions 

Eexp{i&M„[/]} = expl^- l^\9f{xT Mdx)^ . 

In the Gaussian case a — 2, this is just the usual Wiener integral. 

In the rest of this section we develop a discrete approximation of Ma 
when E = with Lebesgue control measure. We begin by recalling 
that the domain of attraction of Sa consists of random variables ^ such 
that 

n 

a~^(^^S_k -bn^ =^ Sa as rn- cx), (5) 
fe=i 

where a„ > and 6„ € M are normalization constants and the S,k 's are 
i.i.d. copies of ^. In the sequel, we will assume r] is SaS , and ^ is not 
only in the domain of attraction of 77, but also that the normalization 
constants are precisely 

On — n^^" and 6„ = 0, n > 1. (6) 

When a < 2, such distributions are said to be in the domain of normal 
attraction of Sa which is not to be confused with the normal domain 
of attraction. 

We propose a discrete approximation of Ala based on the lattice 
hZ"^ C M'^ with edge length h. Let (^fe)fegz<i be a random field of i.i.d. 
copies of ^ satisfying ([5]) and ([6]) and formally define 

where for J'* = [0, 1)'*, : Z'' K is 

f\k) / fix)dx, feLlM')- (8) 

Note that we have implicitly fixed an enumeration {kn,n > 1} of Z'' 
and convergence of X^feeZ'' really means convergence of ^fcn ■ 

The discrete random measures approximate Ma in the follow- 
ing sense: 

Theorem 2.1 (Approximation of SaS random measures). Fix a G 
(0,2]. If a e [1,2], let ft £ L"(M'*) for all t in an index set T. If 
a e (0, 1), for a fixed e > let ft e L""' n Li(M'^) for all t G T. Then 
as h ^ 

M^ift] ^ Maift]. (9) 
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The notation denotes weak convergence of the finite dimen- 
sional distributions, i.e., convergence in distribution of for all 
linear combinations / = + • • • + 9nft„- When the functions are 
indexed by one-dimensional time, it was shown in jAT92| . that even 
for the simple family ft — l[o,t] G the above convergence does 
not hold in the Ji-Skorokhod topologjU. Theorem 1 2. II will follow from 
a Lindeberg-Feller type result for stable distributions which we state 
in Theorem 12.21 below. 

Let us make one more remark before stating Theorem 12.21 One 
motivation for (jH]) was to provide a means to simulate a process Xt = 
Js^d ft{x)Ma{dx). For such simulations, it is natural to let the ^fc's be 
i.i.d. copies of Sa (rather than only in the domain of normal attrac- 
tion). If one is concerned only with one-dimensional distributions (a 
single function /), then a better approximation is given by replacing 
f^{k) in by 

)<l/a> 
(10) 

sign(a::)|a;|". In fact, using 
Uk, one can check that the approximation is exact, and the right and 
left sides of ^ are equal in distribution for every h > 0. The reason 
we have not used (I10|) for the general approximation scheme is due 
to the fact that ([9]) is no longer a SaS random measure under (fTO|) 
because the linearity property ^ does not hold. The analysis of the 
finite-dimensional distributions then becomes much more difficult. 

Theorem 2.2 (Lindcberg- Feller type stable limit theorem). Suppose 
i(.k,j)k,jeN *s an i.i.d. array of random variables in the domain of 
normal attraction of Sa, a G (0,2], and (u'-^')jgN is a sequence of 
vectors in i", i.e. u^'^ := (wi^^)fceN e for all j G N. // 

1. limj^oo II'"'-"'-' life = (7 and 

2. lim^^^oo Wu^^^Wi^ = 

then u'"/^ ^k,j < oo a.s. for each j G N and 



Wfc Cfej =^ <Sa{a-) as j oo. 
ken 



Remarks: 



^ In |AT92j ■ it was also shown tliat under the right conditions, convergence does occur 
in Skorokhod's Mi topology. 
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1. The condition that the S^k.j be identically distributed can be re- 
laxed slightly to the condition that E[exp{i6^k.j)] = 1 — 16*1" + 
o(|0|") holds uniformly in k,j as — >■ 0. For example, they may 
be chosen from a finite family of distributions in the domain of 
normal attraction of Sa ■ 

2. The a.s. convergence J^kefi'^k^ < oo in fact occurs if and 
only if M = (ufe)fegN S as will be seen in Lemma [4. II 

3. Although the series X^feeN^fc '^fcj ^^t converge absolutely, 
switching the order of summation does not change the conver- 
gence in distribution to iSq((t). This will be apparent in the 
proof. 

4. In the Gaussian case, the result can be seen as a variant of 
the usual Lindeberg-Feller Theorem by noticing that condition 
2, concerning loo, is equivalent to 




for all e > 0. More generally when < a < 2, the result is related 
to Theorem 3.3 of [Pct95j which gives necessary and sufficient 
conditions for convergence of sums of independent triangular ar- 
rays to a given infinitely divisible distribution. In particular, the 
conditions of Theorem l2.2l above imply the infinite smallness con- 
dition (cf. Eq. (3.2) in Pet95 ). However, it is unclear how to 
obtain Theorem 12.21 from [Pet95i Thm 3.3] in a manner simpler 
than the proof of Theorem 12.21 provided below. 

3 Linear fractional stable random measures 

To simplify matters, in this section we will restrict our attention to the 
one-dimensional case E = equipped with Lebesgue measure. For 
higher dimensions, see the first remark following Corollarv 13.11 Also, 
in this section we assume that 1 < a < 2. 

3.1 Fractional integro-differentiation and L-FSM in- 
tegrals 

In this subsection we define the stochastic integration of suitable func- 
tions with respect to different L-FSMs in terms of stable random mea- 
sures which are not independently scattered. This is achieved using 
fractional integrals and derivatives. The intuition behind our defini- 
tion is based on two facts. The first is that fractional integrals and 
derivatives can be realized using convolutions, and the second is that 
convolutions are moving averages. 
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The practice of using fractional integro-differentiation for analogous 



integrals with respect to FBM was initiated in |DU99j , and was subse- 
quently used in [PTOO] . We note that the M operator, which is funda- 
mental in the development of the so-called WIS integral f |EVDHd3] ). 
is simply fractional integro-differentiation in disguise. 

Before we define our integral, let us review some preliminaries con- 
cerning fractional integro-differentiation. The Riemann-Liouville inte- 
grals are defined, for / G LP(R), I <p <l/5 and < (5 < 1, by 



^ r fit) 

m J-oo - ty 

1 f fit) 



— / dt (11) 



dt 



ns) Jr ix ~ t) 



l-S 



]_ f fit) 
m Jm ix~t)'_ 



ii-f)i^) — I , .Ls dt (12) 



Our notation is consistent with the standard reference on this topic, 
|SKM87| Sec. 5.1], where some basic properties of the above can be 
found. For example, if / is in the Schwartz space and we allow for 
(5 G N, then PT|) gives the usual integral, as can be seen by Cauchy's 
formula for repeated integration: 

1 



f{h) dh--- dtn^^dtn = -. — / (.T - 0"" V(i) dt. 

{n-\)\]_^ 

Also, the above fractional integrals have the semigroup property for 
(5, 7 > and 5 + 7 < 1: 

44/ = 4+V- 

For sufficiently nice /, this semigroup property extends to all (5,7 > 0. 

Suppose / G and /' G . These are sufficient conditions for 
the following Riemann-Liouville derivatives to exist: 

If / G L^, it is known that the inversion 23£4/ = / holds. 

Bringing the derivative inside the integral in (|13|) . the Riemann- 
Liouville integrals and derivatives of / can be seen as convolutions of 
/ and /' with the family 

WaA^) = y^ali^) ■■= ax-J + hx+^, 13 G (0, 1) (14) 
where we have set (5 — 1 — 5. 
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Definition 2 (Linear fractional stable random measures). Fix 1 < 

a < 2 and a,b > 0. 

1. Ifpe (1/a, 1), let f eL^n L". 

The linear fractional random measure with long range dependence 
is defined by 

M^Mf] ■■= AUaI^_f + b4f]=MM*wi^,h (15) 
where the Hurst parameter is given byH=l + l/a — /3. 

2. If 13 e (0,1/a), let f e and f e flL". 

The linear fractional random measure with anti-persistence is de- 
fined by 

M^Mf] ■■= M^[aVtf + bVlf] = M„[^(/ * w'f})] (16) 

where H — \/a — [3. 

It is not hard to check that I^f and Vj^f are in so that and 
(|16p are well-defined: to see this, split w'^fl into an and function 
using l[-e,e] + l[-e,e]<= ^nd apply Young's convolution inequality, 

Il/*5ll.<ll/llpll5ll, for i + i = i + l. (17) 

p q r 

In fact, one can slightly improve the condition for (fT5|) to / e 7^a(i+a(i-/3)+e) ^ p 
^a(i+a(i-/3)-€) Jqj, some e > 0, and a similar condition can be found 
for (|16p and /'. However, in the interest of simple notation, we will 
not utilize these meager improvements in the sequel. Let us remark 
that the fact that (|15p is well-defined coincides with Proposition 3.2 in 
[PTOO| for the Gaussian case. 

By the linearity of convolutions, it follows that Ma.H is a SaS 
random measure. Also note that Ma^nif] can be interpreted as the 
integral of / with respect to a L-FSM in which case we write 

M,,ff[/] = / fdL^^H- (18) 

To check consistency with we see that 

Mo,[l[0,t] * Wa^b] (19) 

l[o,t](y) {aix - v)-J + b{x - y);^) j M^{dx) 



[o,t](y) (a(y " 2;)/ + b{y - x)J^ dy^ M^{dx) 



ft'\x)M^{dx) 



9 



and 

dx 



M<.[^{\o,t]*Wafi){x)] (20) 



^ l[o,t](y) {a{x - y)J + b{x ~ y)/) dyM^{dx) 

f^'\x)A'Udx). 

When f € and that /' E L^, one can rewrite as 
1 ffi^-t).,, 



r(i - /3) 7k 



+ 



r°° r°° 1 



r(i-/3)7o ^ ^ V, 

- ds. 



r(i-/3)7o 5^+^ 

The right-hand side above is shghtly more general then and is 
called the Marchaud derivative. This is the fractional derivative used 
in jPTOO) ■ however, to keep a unified notation in our approximations 
of the next subsection, we will continue with the Riemann-Liouville 
derivative. 



3.2 Discrete approximations of linear fractional sta- 
ble measures 

Let 1 < a < 2, and consider the stationary moving average process 
(C'c)fcez obtained by "linearly filtering" an i.i.d. sequence in the 

domain of normal attraction of Sa '■ 

a :-5]i;fc-i6- (22) 

Lemma HTTI shows that if u G the series (|22l) converges almost surely. 

Recall the definition of from ([5]) and denote the inversion of a 
sequence by Vk '■= V-k A first stab at approximating a L-FSM integral 
of /, as defined in the previous subsection, might be to mimic (O and 
look at X^fcez fk^k for appropriate filters v (which would also depend 
on h). This is, for example, the approach of |KM88j and |MS08| . Then 
formally, 

J2f^^>^ = E/^'=-'^' (23) 
feez fcjez 
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However, in view of the right-hand side above, it is easier and per- 
haps more natural to first convolve / with Wa,b — w^fl ^'^d then ap- 
proximate the convolution on a lattice with side-length h. In particular, 
for Wa^b corresponding to H E (l/a, 1), define 

:=E(/*^".'')fc^'=' feL'nL"{R) (24) 
feez 

where the sequence {f * Wa,b)^ is defined according to ([S|). Alterna- 
tively, for Wa^b corresponding to H G (0, 1/a), define 

M^mU] --^ E (/' * ^fc' / ^ /' e n L"(M). (25) 

feez 

By ([TT]) and the remark above it, / * Wa.fc £ L"(R). Thus one obtains, 
from a direct application of Theorem 12.11 the following corollary: 

Corollary 3.1. Fix 1 < a < 2. Suppose that for all t in an index set 
T, ft e n L°' when H e (l/a,l) or ft £ C^,// £ n when 
H e (0,1/q:). Then as h ^ 0; 

M^^^ift] ^ M^Mft]- (26) 

Remarks: 

1. It is not hard to extend the H > 1/a case to stable random 
measures on K'* by generalizing the two fixed values a,b > 
(representing the negative and positive directions) to a function 
on the unit sphere S'^^^ C M''. However, one then has to specify 
what is meant by "stationary increments" as there are different 
possibilities for d > 1. 

2. Extending the H < 1/a case to higher dimensions is more dif- 
ficult. One possibility is to consider the Marchaud derivative 
in place of the Riemann-Liouville derivative (see also the next 
remark) . 

3. When H > 1/a, Eq. ([26]) has been shown by various authors 
in the case where ft = l[o,t] (see |KM88j and its references). 
However, when H < 1/a, to our knowledge, even the case ft = 
l[o,t] has not appeared in the literature. It is, however, related 
to the normalization suggested in Theorem 5.2 of |KM88] which 
can be thought of as a discrete Marchaud derivative in the case 
where = l[o ^j. 

4 Proofs 

Before delving into the proofs, let us recall some facts about the domain 
of attraction of a stable distribution. We write f{x) ~ g{x) as x — > c 
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if limx^c f (x) / g{x) ~ 1. For a G (0,2), the fohowing statements are 
equivalent (see |GH971 Theorem 1] with p = 1/2): 

i) ^ is in the domain of attraction of Sa (i.e. Eq. ([S]) holds); 

ii) the tail function t i— >■ P(|^| > t) is regularly varying at infinity 
with index -a and P(^ < -t) - P(^ > t) as t oo; 

iii) the characteristic function X{d) = E [e*^^] satisfies 

- 6 1-^ 1 — Re{X{9)) is regularly varying at with index a, 

- for all X ^ 0, 

a;Im(A(fe)) - Im(A(6l)) 
1 -Rc(A(6i)) ~ 

Moreover, if conditions (i)-(iii) hold, then 

/>oo 

1 - Re(A(6')) - c„P(|C| > l/f) withc„=/ x^"smxdx 

Jo 

as — > and also 

.1/9 

Im(A(6')) = 9 / (P(C > s) - P(^ < -s))ds + o(P(|^| > 1/61)). 



Also, Remark 3 of [GH97) shows that one may choose the normalization 
constants so that 

lim n( 1 — Rc(A(l/a„)) ) ~1 and &„ = nIm(A(l/a„)). 

Recall from ^ that in the present framework, we have assumed 
a„ = n^/" and 6„ — 0. Thus, 

P(C > ~ < -t) - :r~^"" as 6* -j> 0. (27) 

and 

A(6') = 1 - +0(16*1") = A„(6i) +o(|6i|") as 6* ^ 0. (28) 

where Aq. is defined in Eq. Furthermore, (l?7|) implies that there 
exists C > such that for any s > 

Var[a{|c|<.}] < and E[|e|l{ici<«}] < Cs^-". (29) 

4.1 Proof of Theorem [272] 

We begin with a lemma which shows the is the right space for the 
sequence u. 
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Lemma 4.1. If {^k)kefi '^'^ i.i.d. sequence in the domain of normal 
attraction of Sa, then X^fe'^^feCfc < if and only if u ^ . 

Proof of Lemma \4-l\ The case a = 2 is standard and omitted. Con- 
sider a € (0,2). Recall Kolmogorov's Three-series Theorem: J2k''^k£,k 
converges a.s. if and only if for any s > 0, the following three series 
converge 

fceN fcsN fceN 

Eq. (gT)) implies 

and hence the first series converges if and only if u S If m e i?", 
then ([^5]) implies the convergence of the third series since 



\uk\mk\h\e,\<s/M}] < \uk\c{s/\uk\r-' 

The convergence of the second series is obvious. 



= Cs 



Uk 



□ 



Proof of Theorem \2.2[ If rjkj are i.i.d. iSq random variables, then for 
any fixed j the above lemma allows us to write 



Eexp{i9Y,4'^Vk,j} ^ II K {u''^^9^ =exp{-|||u(-')l|,„0p} (30) 
and 



ken 



keN 



Eexp{z0^ = n A (^^^fl) . (31) 

fceN fceN 

Note that since := J2keN !'"(''') I" absolutely converges, the order 

in which the summation and products above are taken is irrelevant. 
It suffices to show that as j oo. 



(32) 



keN 



ken 



We fix j and estimate the difference of the above products using the 
following fact : if {zi)i^i and (zi)ie/ are two families of complex num- 
bers with moduli no greater than 1 and such that the products Yiiei ■^^ 
and Hie/ -^i converge, then 



<^\zl-Z^\ 



iei 



iei iei 
We therefore have 



(33) 



fcGN 



kef 



ker 



A, (ui^h) I . (34) 
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Equation (|28)) impliefQ that the function g defined by 17(0) = and 

g{u) ^ \u\-" \X{u) ~ X^{u)\ , u^O, 
is continuous and bounded and for any A: G N, we have 

In order to obtain a uniform estimate on the above, define the function 

g : M+ ^ R+ by 

g{v) := sup \g{u)\. 

\u\<v 

Note that g is continuous, bounded and vanishes at 0, and that for any 
fc e N such that \u[^^9\ < £, 

x[u\^^9)~X^ [u^^^e)\<~gie)\u^^^9r. 

Let £ > 0. Equations and psp together yield 



(35) 



fceF 



kef 



fcGF 



9|>£}- 



Now, by the continuity of g at 0, g{e) is smah when e is smah. Eq. 

follows since lim^^oc = implies I]fceN l{|„0'e|>e} ^ ^ 

as j — 00. □ 

4.2 Proof of Theorem 

Let [-J denote the floor function applied to each coordinate of W^. 
Define fh : M'' t-^ M to be a piece- wise constant function approximating 

/ e a 



loc\ 



h{x) 



h-'fiy) dy 

h([/l-lIJ+/<*) 

h'^^fiv) dy, for xeh{k + I'^) 
for X e h{k + I'^). 



(36) 



h(k+I<i) 



Note that 



(37) 



'We have assume a £ (0, 2) for Eq. (|28|) . but for o = 2 it is well-known. 
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Lemma 4.2. For a G [1, 2], suppose f E L"(]R''). Then as h ^ 0, 

lim =0. 

Proof of Lemma \4.S\ To reduce notation we assume d = 1, but the 
proof holds for general d. Fix fc 6 Z and consider the sequence of /I's 
such that h — 2^^ for j e N. We will exploit the fact that fhi[k,k+i) is 
a martingale (in time j) with respect to Lebesgue measure on [fc, fc + 1) 
and with respect to the cr-fields generated by the sets i + i G Z. 

For a > 1, |/h|"l[fe,fc+i) is a submartingale which, by the martingale 
convergence theorem, converges a.s. to |/|"l[fe,fe+i) ■ Thus, Fatou's 
lemma gives 

||Al[fc,fc+l)||L° ll/l[fc,fe+i)llL°- (38) 

Since fh^[k.k+i) converges a.s. and the L"-norms converge, we have 
convergence in L"(R) of //il[fc,fc+i) and also for //il[_Ar.Ar) for any N g 
N. 

For / £ i"(M) without compact support, simply choose N so that 

Wf^-N^NyWl" < 

Since |//i|"l[fc,fc+i) is a submartingale, we also have ||//!.l[-Ar.Ar)c < e 
uniformly in h. 

Finally, to extend the above to general h 0. Note that all we 
really require is a sequence of lattices such that finer lattices are sub- 
lattices of prior ones and that the mesh size goes to zero. But any such 
sequence has the same limit in L"(R), thus we conclude that the only 
real requirement is that the mesh size goes to zero. □ 

Proof of Theorem \2.1i By the Cramer- Wold device, we must show that 
for aU 01, 0„ e M and /i, ...,/„ e L"(K'^), 

n n 

OrM^ [/,] =^ J2 as /i 0. 

i=l 1=1 

Our proof uses Theorem 12.21 First note that the comment following 
PTjl shows that switching the order of summation in the series [fi] 
does not affect its distribution. This, together with the linearity of Ma 
and M^, allows us to reduce the above to verifying 

M^[f]=^MM] ash^O 

for a sing le / e L"(R''). This will follow from Theorem O provided 
we check the two conditions 

hm ll/'^ll,. = ||/|U» (39) 
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and 



lim 11/^1 



0. 



(40) 



We consider a E [1,2] first. Condition easily follows from ([57)) 
and Lemma [4.21 For (|40p . note that convergence of the norms of 
coupled with a.e. convergence, shows that the family {\fh\"}hGN 
is uniformly integrable. 

For a G (0, 1), we first consider the sequence of /I's such that h = 

for j € N. By uniform integrability and the martingale convergence 
theorem (see the proof of Lemma . we see that fh^[k.k+i} converges 
in to fl[k,k+i)- The final comment in the proof of Lemma 14. 2 1 shows 
the convergence also holds for arbitrary h ^ 0. 

Next, note that L^{[k, k + 1)) contains L"([fc, k + 1)) and that the 
endomorphism on {[k , k + 1)) which maps 

/1[/C,/C+1) ^ |/ri[/£,/£+l) 



is continuous. Thus Eq. (1551) holds for a G (0, 1). 

Since / S we can choose Ni so that ||/l[-Ari,Ari)<: ||lc is small. 
However, to uniformly bound the tails of the fh, we will use the stronger 
condition of / G L"""^. In particular, there exist N2 > 0, C > and 
S > such that \x\ > N2 implies |/(a:)| < C\x\~^ . We have for 
|a;| > N2 + h that 



and 



h-' I f{y)dy 

h(lh-^x\+I} 



< Ch^-^ilxl-h) 



-aS 



(41) 



Since aS > 1, ([M|) follows from ([55)) . Finally, as before, we see that 
([5^ along with a.e. convergence gives ([^ for a S (0, 1). □ 
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